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We study the Directed Polymer model subject to a particular form of disorder, rj(x, t) — rix(x)r]T(t), 
recently proposed in biological applications. We find that two new universality classes arise, depend- 
ing on the the lattice geometry. Using an intermediate model linking the two different orientations 
continuously, we find that there is a phase transition separating two distinct scaling phases. For both 
phases we get a reasonable understanding of the nature and values of the exponents, corroborated 
with numerical results. 
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In their study of DNA-sequence alignments, Hwa and 
Lassig consider the mapping of the problem of sequence 
recognition and matching to a more standard (in statisti- 
cal mechanics) model of directed polymers (DP) subject 
to a special disorder Q . They establish, both analytically 
and numerically, the equivalence of the two problems. In- 
deed, their model is extremely similar to the DP problem 
H but with a special disorder structure. Whereas in the 
usual DP problem every site (x, t) on the lattice has an 
energy 77(01, t) independent from the energy of any other 
site, Hwa and Lassig define 77(2, t) as 

f](x,t) = 7] X (x)r] T (t) , (1) 

where {r]x(x)} and {?7t(£)} are independent unbiased 
(i.e. with zero average) random variables. For the bi- 
ological motivations to consider this kind of disorder we 
refer the reader to Hwa and Lassig's work [Q . Naively one 
may expect that this disorder distribution should give no 
difference with respect to the usual uncorrelated disorder. 
Indeed correlations may appear non-relevant, since 

77(3:1, h)r)(x2, t 2 ) = r)x{x 1 )f}x{x^) ifr{ti)T]T(t2) (2) 

is zero if x± 7^ X2 and t\ ^ as for the usual uncorre- 
lated disorder. However, the count of the total number 
of disorder degrees of freedom suggests that the model 
should be very different. On a 2-dimensional square lat- 
tice of linear size L (L 2 sites) there are only 2L degrees 
of freedom, versus L 2 degrees of freedom for the usual 
site disorder. In other contexts this difference leads to 
relevant and important consequences. Indeed, we recall 
the Mattis model of spin-glasses 0j , where spin coupling 
disorder is replaced by a product of disorders, and frus- 
tration effects disappear. 

In this Letter we systematically study the effects of 
such a factorized disorder on the properties of directed 
polymers, both with numerical simulations and with the- 
oretical arguments. We find that two new universality 
classes emerge, depending on the orientation of the lat- 
tice. We then show that it is possible to obtain both 
universality classes on a same system by varying a pa- 
rameter, a signature of a phase transition. Finally, we 



give explanations for the difference between our results 
and Hwa and Lassig's ones 0] . 

Eq.(Q) defines the model on a rectangular lattice with 
disorder defined along the (1,0) and (0,1) directions. 
We can also define the model on the diagonal (1,1) and 
(—1,1) directions, i.e. using the so-called diamond lat- 
tice. In the latter case Eq.(jl]) becomes 

r)(x,t) — 77(7; + x)i](t — x) (3) 

as was originally studied in jj] . 

First we study numerically the model on a diamond 
lattice, eq.(|J). As for the usual study of the ground 
state of directed polymers (i.e. temperature T = 0), 
we are interested in the scaling properties of the poly- 
mers as a function of their length. In particular we con- 
sider the transverse wandering fluctuations of the poly- 
mer, < \5x\ >~ t v , where angular brackets indicate the 
average over different realizations of disorder (as usual 
in quenched disorder problems). Also, the ground state 
energy fluctuations are of interest, < \SE\ >~ t u . 

We always use disorder taken from uniform distribu- 
tions over compact supports, in general in the interval 
[—1/2, 1/2]. We discuss later the effects of biased distri- 
butions. 

A summary of the exponents obtained is given in Table 

I 

For the d = 1 + 1 diamond lattice we observe indeed 
the scaling of v = 2/3, consistent with the standard dis- 
order result (|], see Figj^. Yet, the ground state energy 
fluctuations scale with an exponent u> = 1/2 different 
from the standard w = 1/3. In order to appreciate the 
full mechanism of the model, and to understand where 
the difference comes from, we study the time evolution 
of the energy profiles on diamond lattice strips of various 
width L. It corresponds to the growth of a KPZ-like 
surface in discrete space. 

The beahvior of the energy profile fluctuations 

AE{L,t) =< (E(x,t) - E(x,t)) 2 > x ' 2 ( over-bars rep- 
resent averages over one single disorder realization) is 
governed by three exponents. Initially AE(L, t) grows 
with time as AE(L,t) ~ t@. Then it saturates to a con- 
stant value, dependent on the size L of the system as 
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AE(L,t) ~ L x . The characteristic time r separating 
the two regimes scales with the size L of the system as 
r ~ L z , and z is called dynamical exponent. Scaling con- 
sistency imposes z = x//3. In standard DP-KPZ theory 
the surface exponents and the polymer ones are related 
by z = 1/v and (3 — u>. Moreover, in standard KPZ the- 
ory the law z + x = 2 holds, due to galilean invariance. 

The scaling analysis of our simulations gives \ = 1/2, 
as in standard KPZ theory, and (3 = 1/2. From scaling 
consistency we have z = x//3 = 1, which is in agree- 
ment with numerical simulations, but that does not sat- 
isfy neither the relation z = 1/v = 3/2, nor z + \ = 2. 
The reason behind the above scaling behavior, namely 
z = 1 instead of z = 1/v = 3/2, is due to the particular 
structure of the disorder landscape. Indeed, the diag- 
onal (±1,1) factorization of disorder imposes a space- 
time linear proportionality relation, corresponding to a 
dynamical exponent z = 1. We therefore find that the 
usual relation z = 1/v breaks down. Let us now try to 
explain the 1/2 value of the x an d u> = (3 exponents. 
Indeed the u> = 1/2 could come from an uncorrelated 
random walk, whereas x = 1/2 is consistent with the 
standard result. In order to solve this ambiguity we have 
analyzed the energy profile evolution for a 2 + 1 dimen- 
sional diamond lattice. Again we find 2 = 1, strength- 
ening the time-space proportionality interpretation, and 
(3 = x = 0.40 ± 0.01, nicely consistent with the standard 
result We conclude therefore that the asymptotic 

behavior is KPZ dominated. At variance with standard 
KPZ, the approach to the stationary state changes, gov- 
erned by a ballistic dynamical exponent z = 1. 

Let us now consider the square lattice with disorder de- 
fined from Eq. (|lj). The directed polymer scaling results 
are shown in Fig.^|. The exponents are consistent with 
v = 2/3, as in the diamond case, and to = 2/3. From 
energy profile evolution, we find the exponents (3 = 2/3 
and z = 3/2 (consistent respectively with (3 = lo and 
z = 1/v)- The roughness exponent is then x = z/3 = 1, 
in agreement with numerical results. The absence of any 
space-time relation imposed by the disorder structure 
leaves therefore the usual z = 1/v relation unchanged. 

We explain the value 2/3 of the v and lo exponents 
using an argument generalized from columnar disorder 
theory HQ . Indeed, as emerges from FigJ^, after an ini- 
tial wandering, polymers tend to localize in a very nar- 
row region of the lattice, in general consisting of a pair 
of columns. The corresponding two r\x(x) values have 
opposite signs and their absolute values are extremely 
close to 1/2. Polymers are then attracted to such pairs, 
because they can jump from one column to the other 
depending on the sign of ifr(t), always choosing the neg- 
ative product rjx{x)rjT{t) in their quest for low energies. 
We can therefore think that polymers search the best 
pairs available. The relevant feature of the system is 
therefore the effective pair energy probability distribu- 
tion Q(e). From its knowledge, using variational argu- 



ments, it is possible to obtain an analytical estimate of 
the exponents. Indeed, since the polymer always chooses 
the negative product, as long as we are interested in an 
effective pair energy distribution we can disregard the 
rjT{t) energies, just remembering that the chosen rjx(x) 
energies are extremal (i.e., as close as possible to ±1/2). 
Moreover a non-restrictive extra condition is that the two 
r]x{x) energies are equal in absolute value. This is nu- 
merically well approximated, since the two energies are 
both extremely close, in absolute value, to 1/2. 

Let us define e such that the energies of the two rjx (x) 
of the pair are ±(1/2 - e). Then e e [0, 1]. The prob- 
ability that the pair energies are closer than e to ±1/2 
is 



(4) 



leading to a probability distribution for e that is Q(e) = 
2e. Q(e) is then the relevant pair effective distribution. 

It is then possible to apply a variational argument to 
give a prediction on the exponent v. Indeed, the polymer 
wanders initially in a region of characteristic width R, in 
search of the optimal (pair) energy. Once it has found 
the optimal pair, it localizes there up to a time t. The 
global energy can therefore be evaluated as 
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E = 1 R+(t- R)e r 



(5) 



The optimal energy available in a region of width R can 
be obtained from the statistics of the extrema. Indeed, 
the probability Q m (e) of the smaller energy among R 
energies distributed according to Q(e) — (/i + l)e' i in the 
interval [0, 1] is 



i- / Q(y)dy 
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(6) 
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Q m (e) = RQ(e) 

i?( M +iy (l-e^ 1 ) 

Then the characteristic extremal energy e m i n is 

e min = J Q eQ m (e)de = Rb(r- , (8) 

where B(x, y) is the Euler Beta function. After some 
manipulation, the scaling behavior of e m i n is found to be 
£min ~ for large R. Inserting this result in 

(||) and optimizing with respect to R, we find the scaling 
behavior R - i^ 1 )/^ 2 ) an d therefore v = 
2); in our case /i = 1 and therefore v = 2/3, as from 
numerical simulations. Also, the energy E scales as R 
and therefore uj = v, again in agreement with numerical 
simulations. 

We can also change the probability distribution of the 
energies along the x and t directions. As long as both en- 
ergies can take on positive and negative values, and the 
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polymer does not loose any energy wandering through 
different columns, the universality class does not change. 
It does not even change if the energies ijx (x) take only 
positive values: the polymers continue to localize on pairs 
of energies, to maximize the energy gain when the r]T{t) 
energy is negative and to minimize the energy loss when 
r)T(t) is positive. The only different universality class 
comes when the rfr(t) energies assume only positive val- 
ues: in that case there is no more any convenience for the 
polymers to bind to a pair, and localization takes place 
over a single optimal rjx {x) energy. Being the energy dis- 
tribution flat (/x = 0), the exponents are v = u = 1/2, in 
agreement with numerical results. If instead there is an 
energy cost to change column, then a transition as a func- 
tion of the bias in the distributions takes place between 
the two scaling regimes v = 1/2 and v = 2/3. The full 
phase diagram of the latter situation will be published 
elsewhere. 

This polymer problem defined by (|l]) is indeed sensi- 
tive to the chosen disorder distribution. For instance, 
using Gaussians distributions, the variational procedure 
predicts a linear relation between R and t (apart from 
logarithmic corrections fft), confirmed by our numerical 
simulations. 

Next, we move on to show that the two different scaling 
behaviors (namely, z — 1 and z — 3/2) are compatible 
on the same disorder configuration. Let us start from 
the square lattice: a polymer has three possibilities to 
go one step forward, namely, to the left, in front, to the 
right. In our simulations, we chose total isotropy, that is, 
neither of the three choices imply any energy loss, and 
the z = 3/2 result is recovered. We let then one of the 
directions, say the leftward one, be energetically unfavor- 
able, with an associated energy cost a. Consequently the 
polymer end-points show a drift in time toward the right 
direction. Indeed, we find that, for very large values of 
a, the polymer average final position moves to the right 
linearly with time. The fluctuations of the end-points 
around their average position scale with time with expo- 
nent v = 2/3, and the energy fluctuations scale with ex- 
ponent u) = 1/2, as for the diamond lattice case. Indeed, 
looking at the scaling properties of the energy profile, 
we recover the full set of exponents typical of polymers 
with factorized disorder on the diamond lattice, Eq.([|). 
Strong anisotropy induces therefore a linear space-time 
relation (as emerges from the linear drift of the polymer 
end-points) that forces a z = 1 dynamical exponent. For 
small values of a (down to the isotropic case a = 0), in- 
stead, we find the other full set of exponents, typical of 
disorder (1). 

As Fig. 4 shows, by varying a we find the signature of 
a phase transition between the two regimes z = 3/2 and 
z = 1. The three point sets are obtained from collapse 
plots of the scaling features of energy surfaces for different 
system sizes. Precisely, the values plotted come from the 
collapse of data of couples of strips. The values of the 



width pairs are given in the legend. The data are steeper 
with increasing strip width in the neighbor of a = 0.25. 
We expect a = 0.25 to be precisely the transition value of 
the anisotropic energy. Indeed, localization takes place 
when, on the average, it is convenient for polymers to 
localize even if it has to pay a price a: a localized polymer 
chooses the column of the two where it finds negative 
disorder energies. Since the sign of rjx(x) is fixed, it finds 
a sequence of n negative values on the same column if 
r]T(t) keeps the same sign n times, and this happens with 
probability (1/2)™. The average length n of a sequence 
with the same sign is therefore 



71=1 



(9) 



The average energy \r)T{t)\ = 1/4 (the sign depends on 
the sign of rjx(x)). The absolute value of the energy 
r]x{x) is extremely close to 1/2 (again irrespective of the 
sign, which depends on the sign of ryr(t) so that their 
product is negative). The average energy that a localized 
polymer gains against a is therefore 
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(10) 



Therefore if a < 1/4 localization takes place. Such a 
transition was also alluded to in Ref. || . 

Summarizing, we have found that a factorized disorder 
distribution such as (|l|) or (H) leads to new universality 
classes, depending on the orientation of the lattice. We 
have been able to explain the origin of the characteristic 
exponents of these classes and to connect them to each 
other via an anisotropy driven phase transition. 

A challenging issue is to connect these results to Hwa 
et al. [fj]j^], where the usual disorder universality class 
was obtained Indeed, in their work they defined the 
disorder and the lattice in a slightly different fashion (see 
Refs. |l],§] for details). As we showed above, the polymer 
problem with factorized disorder is very sensitive to the 
details of the disorder and of the lattice. We believe that 
the reason underlying this difference resides in the pres- 
ence, in their work, of some directions along which poly- 
mers do not feel any disorder (the gaps in the language 
of Refs. |i],D)- Indeed, depending on the energy cost of 
these directions, we found indications of a phase transi- 
tion between the standard disorder universality class and 
the diamond lattice one (see Table Q). Further work is 
needed in this direction, to explore the full phase diagram 
of the model. 
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FIG. 1. Wandering and Energy Fluctuations for polymers 
on the diamond lattice; data are obtained after averaging over 
10000 disorder realizations. 



1000 



100 



10 





.> 


10 100 1000 10000 A ' 

r -"V »-* ^ 
y » - ■ s 

y W • -"" 




♦ ♦ Wandering 

• • Energy 



10 



100 
Time t 



1000 



10000 



FIG. 2. Wandering and Energy Fluctuations for polymers 
on the diamond lattice; data are obtained after averaging over 
10000 disorder realizations. Roughness of the energy profile 
on strips of different width L, again averages are over 10000 
samples. 




FIG. 3. Ground state polymers on a square lattice. Every 
polymer corresponds to a different disorder realization. 




FIG. 4. Exponent z vs. anisotropy energy a, as obtained 
from collapse plots of roughnesses for different strip widths. 
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UJ=f3 


1/2 


2/3 


1/3 


z 


1 


3/2 


3/2 


X 


1/2 


1 


1/2 



TABLE I. Characteristic exponents for directed polymers 
and energy surfaces for the diamond lattice (disorder defined 
from Eq.(3)), square lattice (Eq.(l)), and for the standard 
disorder [3]. 
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